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A new non-perturbative framework for many-body correlated systems is formulated by extend- 
ing the operator projection method (OPM). This method offers a systematic expansion which 
enables us to project into the low-energy structure after extracting the higher-energy hierar- 
chy. This method also opens a way to systematically take into account the effects of collective 
excitations. The Mott-Hubbard metal-insulator transition in the Hubbard model is studied by 
means of this projection beyond the second order by taking into account magnetic and charge 
fluctuations in the presence of the high-energy Mott-Hubbard structure. At half filling, the 
Mott-Hubbard gap is correctly reproduced between the separated two bands. Near half filling, 
strongly renormalized low-energy single-particle excitations coexisting with the Mott-Hubbard 
bands are shown to appear. The significance of the momentum-dependent self-energy in the 
results is stressed. 
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In spite of a long history, a satisfactory theoretical description of strongly correlated electrons 

remains open because of the difficulty in reproducing their rich hierarchy structures in energy and momentum. 
The properties near the metal-insulator transitions (MITs) provide a typical challenge. It is required to 
correctly describe both a high-energy incoherent structure of the "Hubbard bands" with the Mott gap 
and the low-energy coherent excitations. The low-energy part must be severely renormalized with strong 
momentum dependence. It is also required to-treat magnetic and charge fluctuations. 

The dynamical mean field theory (DMFT)iJlP reproduces both_thc incoherent Hubbard bands and the 
quasiparticle structure. On the other hand, numerical simulationsEr have revealed an important aspect of 
the MIT in the Iwo-dimensional Hubbard model, i.e., nearly dispersionless fermions around the (n, 0) and 
(0, 7r) moment a. 0U This suggests that the momentum dependence of the single-particle self-energy, which 
is completely neglected in the DMFT, is significant in real finite-dimensional systems. The momentum 
dependence of the low-energy structure appears from collective modes under the influence of high-energy 
incoherence. Such hierarchy can not be properly treated in the perturbative or one-loop analyses of the 
self-energy and the dynamical susceptibilities, because the high-energy structure such as the Mott-Hubbard 
bands is missing. 

In this letter, we propose a systematic expansion which incorporates both the coherent and incoherent 
structure in the energy hierarchy of single-particle excitations. By applying the operator projection method 
(OPM) to the single-particle Green's function, our framework allows us to systematically project out the 
high-energy structure and sequentially extract the low-energy coherent excitations. The low-energy coherence 
appears with the strongly momentum dependent selfienergy due teutwo-particle collective fluctuations. _ . 

The OPM first developed by Nakajima,EP ZwanzigeiJ and MoritiP is based on the equation of motiorO 
with a moment expansion of a correlation function. The main procedure of this method is to decompose 
the time-derivative of an operator into the term projected onto itself and the remaining new operator. 
The OPM involves repeatedly applying this procedure to the newly created operators. Then, the OPM 
allows one to construct a continued-fraction expansion of a correlation function. The expansion is suited 
for systematically evaluating the energy hierarchy structure from high to low energies. Similar equation- 
of- motion (EQM) approaches supplemented with various decoupling schemes a n d /or .at h pr a pprmH m a t i nnq 
have bppji intensively studied for strongly correlated models such as the HubbardaEjoli^ and the t-J 
modelsEj'EZlB For the Hubbard model, the Hubbard bands are reproduced from the first three, i.e., the 
zeroth-, first- and second-order moments of the Green's function.© In principle, the low-energy band appears 
in the fourth order, although a systematic and correct expansion has not yet been carried out. In this letter, 
we show that the OPM up to the second order, supplemented with the self-energy in the second-order 



projection which is referred to as the second-order self-energy, makes it possible to correctly reproduce the 
hierarchy: the Hubbard bands as well as the renormalized low-energy bands. We note that our calculated 
second-order self-energy gives an evaluation of the fourth-order moment. The OPM based on the equation 
of motion is a non-perturbative treatment and works from the weakly correlated limit U(n)(2 — (n))/t — » 
to the strongly correlated limit U (n)(2 — (n))/i — > oo. 

We take the Hubbard Hamiltonian with an electron transfer ts.s> from an atomic site x! to x and the local 
Coulomb repulsion U; H = - J2s,x',s ts,S'c' Ss cg' a + U J2s n ^ n n witn n xs = Cg s C£s, or H = H - (i^f s n xs 
with the chemical potential fi. We perform the projection procedure for the electron creation and annihilation 
operators at an atomic site x with a spin index s, ct and cg >s . It is useful to define the fermionic 4N- 
component vector operator composed of 
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Here, N is the total number of the atomic sites. is defined as the Hermitian conjugate operator of \&. 
Co (Po) and o (p) are the identity and the Pauli matrices, respectively, which operate to the spin (charge) 
space of vector operators as ^> from the left-hand side and those as ^ from the right-hand side. 

We introduce CjA = [A,H]-, and (A) = Tr[e~ H / T A]/ r Ti[e~ H ^ T ] at temperature T. We also define the 
response function of *S? in the 4N x AN matrix representation K^^(t). Its (a, a') <S> {x,x?) component is 
defined by 
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We note K^, ^ (0) = —il with the AN x AN identity matrix J because of the fermion anticommutation 
relation. The Fourier transformation of Eq. (|^) gives the susceptibility *t(^), which is nothing but the 
Green's function G(u>). The projection procedure is defined as 



,(n) 



(3a) 
(3b) 



Namely, the operator cD\& is split into with a component given by the equal-time correlation KQWt{0) 
and the new operator Sd)^f which satisfies K^^gQ^y (0) = by definition. It allows us to express the 
single-particle Green's function G(u>) — x<j/ *t(<^) i n the form of a Dyson equation as 
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G(u) = i 
G^M = 
S 1 ( W ) = -ixS*,(i«a*)tMif;^t(0). 
The irreducible part of X<S(i*.(5£*)t ( w ) with respect to G^°\uj) has been introduced; 

XS*,(»4)t( w ) 

= x^^t-V) - ^;,*t(o)G (0) Hl _1 . 



For the Hubbard model, one obtains 
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For the repulsively 
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^(cs|Cjx - ct^ct T ) The summation over 
interacting model (J7 > 0), we can exclude the superconducting states with the isotropic s-wave pairing 
symmetry, (Al^ 4 ) = 0. The thermal average (SW), if it does not vanish, is chosen to break the symmetry in 
the z-axis, i.e., (Si) = 0. Then, Eq. (§) and to* are reduced to 
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5u>^ s = U ( -Sn s a p 3 - 6S~a 3 p 



(9b) 

with Sn s = ng- (rig) and 5S§ = £§- 

If £i(ui) is ignored, then the present theory is reduced to a Hartree-Fock theory. In order to calculate 
S±(lo), one can adopt perturbation expansions on a single-particle basis, one-loop approximation, long- 



time or short-time approximations. However, if Si(uj) is self-consistent 
frameworks such as_Derturbation expansion or one-loop approximations 



mined by means of simple 
W the MIT at half filling- 



can not be reached.EJ This_is because arguments along this line always fail to give the correct high-energy 
behavior of the self-energy r3 We note that only the zeroth- and first-order moments of the Green's function 
are correct in these theories. As known in the Hubbard and two-pole approximations, we need to proceed 
to higher-order moments, as is discussed below. 
The second-order projection gives 
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From Eq. @ one obtains the irreducible self-energj 
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For the repulsive Hubbard model, 
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where // 2 ' = 
^Tre^^o + /O3) : 



U(l — ^(ng)). The kinetic energy per site summed over spins is expressed as £7km — 
— jj s % c fcs)' Below, a momentum independent energy shift £ cor2 due to two-site 



correlated hopping terms which take the form t s g(ct s cg s ,Sng- s ) in Eq. (12b) are neglected. Then, in the 
case of (Sg) — 0, the remaining operator becomes 
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Even when, X l 2 (w) is neglected, the present formalism is a generalization beyond the Hubbard I ap- 
proximations'—because not only the Mott-Hubbard bands are reproduced but also e' 22 ) depends on the 
momentum.liijl!-l ) The former fact is guaranteed, since the Pauli exclusion principle is satisfied. The latter 
fact modifies the single-particle dispersion in each Mott-Hubbard band. The momentum dependence of e^ 22 ) 
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enters the present theory through the equal-time charge, spin, isotropic s-wave pairing correlations, while 
these correlations are ignored in the Hubbard I approximationEP In the symmetry-unbroken phase, the 
electron Green's function with momentum k and spin s reads 
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Here, is a (1,1) component of Eq. (12c). Equation ([14]) has a pole 
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£g) 2 is the momentum-dependent Mott-Hubbard gap. The momentum 
dependence diminishes as U ft diverges. It is easily found that when the collective degrees of freedom are 
completely frozen, i.e., = 0, Eq. (|l4|) recovers the Hubbard I solution at half filling. On the other hand, 
when equal-time two-particle correlations with the momentum Q are significant, behaves like t^+Q' Then, 
Eq. ([l4]) smoothly connects the symmetry-broken and symmetry-unbroken phases through the growth of 
short-range correlations in Eq. (12c). For example, it reproduces the shadow bands in systems with well- 
developed short-range antiferromagnetic correlations. It can also describe other kinds of shadow bands 
such as nearly ferromagnetic systems. So far, few theories iaite reproduced such shadow-band effect even 
qualitatively within the one-loop approximations to X'i(a;).Ej'LfP Without 27a(w), the expectation value of 
the Hamiltonian as a function of is given by 
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with a fermionic Matsubara frequency uj n = (2n + 1)ttT and the Fermi distribution function f(u>) = (e UJ ^ T + 
1) . The momentum-fc spin-s component of £\{ui) has been introduced as Yli e>s (u>,k). At (n) = 1, we 
obtain (H) — — $^r[<£ — (*£ + tj:) 2 /4]/U when U ^> T, t. This includes the superexchange term through Eq. 
(12c). ^ k 

Next, we consider the second-order self-energy £2{u) by introducing decoupling approximations. Here, we 
can take several methods to calculate .^(u;): (A) the perturbation theory, IB) one- loop approximation with 
two-particle susceptibilities obtained by the self-consistent RPA (SCRPA),E3 and-XC) one-loop approxima- 
tion with those obtained from the two-particle self-consistent (TPSC) method. E3'E-iP In the TPSC method, 
irreducible vertices are determined so that the local-moment sum rule is also fulfilled for th£_two-particle 
properties, not only charge and spin correlations^ but also local s-wave pairing correlations.cS Therefore, 
we take the last method, (Cl_We note that ^(w) has been evaluated within a two-site level by means of a 
different decoupling scheme .113' However, short-ranged correlations are ignored there. 

Then, the second-order electronic self-energy is evaluated by means of the decoupling approximation as 
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2nmT is a bosonic Matsubara frequency, and % c , \s and 



Xp are the charge, spin, local s-wave pairing susceptibilities, respectively. The spin indices have not been 
included since all the variables are independent of the spin index. Within the TPSC approximation, the 
susceptibilities are calculated from 



Xc(^m,fc) = 2x P h(f2m, k)/[l 
Xs(tt m ,k) = 2x P h(^m, k)/[l 
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Xp(^m, k) =Xpp(^m, k)/[l + r p Xpp(^m, k)], 



(16c) 



with an approximation T s = AU (n£tn%\) / (n) and the self-consistency conditions for T s , T c and T p given by 
the local moment sum rules; 



T / N ^2mX^( n "i,k) = (n) +2(n s ^n si ) - (n) 2 , 
T / N ^2 m ,kXs(tt m ,k) = (n) -2(n s1 n sl ), 
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Xph and x PP are the particle-hole and particle-particle susceptibilities calculated from the bare Green's 
functions. 

The calculated density of states for the Hubbard model with U = it and the electron transfer restricted to 




Fig. 1. The local density of states for U = At and T = 0.024 at n = 1 (solid line) and n = 0.9 (dashed line). 




Fig. 2. Single-particle spectral weights A(u>, k) = -ImG C|S (w, k)/ir for U = it and T = 0.02t at n = 0.9. 



the nearest neighbor at the fillings n = 0.9 and n = 1 are shown in Fig. 1. The calculations have been done 
in the 32 x 32 lattices with 2048 Matsubara frequencies. The metal and the Mott insulator are reproduced 
at n = 0.9 and n = 1, respectively, with the Hubbard bands. The emergence of a flat dispersion around 
(tt, 0) is in agreement with the numerical results.tf It is remarkable that at n = 0.9, low-energy peaks around 
the magnetic Brillouin zone lie below the chemical potential and the dispersion outside the zone hecomes 
strongly small and pinned near lu = 0, as shown in Fig. 2, which agrees with previous results E'H It 
turns out that the growth of the equal-time antiferromagnetic (AFM) spin correlation produces such strong 
modification of the dispersion in the scale of t. At n = 1, t/t should become —1 at the zero temperature 
to yield the antiferromagnetism with the broken translational symmetry In the present calculation carried 
out at T = 0.02<, in fact, t/t = -0.84 for n = 1, while t/t = -0.53 for n = 0.9. This indicates that as 
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holes are doped into the Mott insulator, the shadow-band dispersion moves toward the chemical potential. 
However, self-consistent calculations of -£2(0;) using the TPSC method tend to weaken the SDW shadow- 
band structure too much and to overestimate the low-energy flat bands, as shown in Fig. 2, because the 
TPSC method overestimates S2 e (w, k) at low frequencies. This will be improved by considering higher-order 
projections, as is discussed elsewhere. 

We also note that in the results calculated from the (A) or (B) method, the AFM spin correlations appear 
to be underestimated. These give qualitatively similar results but smaller values of \t\. 

In summary, the OPM has been applied to the two-dimensional Hubbard model. Beyond the second=order 
projection, which correctly reproduces the first four moments and yields the two-pole approximation,!!^ the 
higher-order dynamics in the self-energy have been self-consistently taken into account. Then, we obtained 
not only the Mott insulator at half filling but also nearly dispersionless low-energy excitations pinned near 
the Fermi level away from half filling. This is the first analytic theory which has succeeded in describing 
both of the Mott-Hubbard bands and the low-energy single-particle excitations near half filling under strong 
short-ranged spin correlations. The OPM applied to the two-particle susceptibilities beyond the SCRPAEj' 
are now under study to give more correct two-particle properties. Higher-order projections are necessary_to 
consider the non-local pairing correlations, which are crucial in relation to the high-T c superconductivity .c^P 
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